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Abstract 

The caustics of Fourier integral operators are defined as caustics of the corresponding Schwartz kernels 
(Lagrangian distributions on X x Y). The caustic set E(C) of the canonical relation C is characterized 
as the set of points where the rank of the projection tt : C — > X x Y is smaller than its maximal value, 
dim(X x Y) — 1. We derive the L P (Y) — > L q (X) estimates on Fourier integral operators with caustics 
of corank 1 (such as caustics of type A m +\, m € N). For the values of p and q outside of certain 
neighborhood of the line of duality, q = p' , the L p — » L q estimates are proved to be caustics-insensitive. 
We apply our results to the analysis of the blow-up of the estimates on the half-wave operator just before 
the geodesic flow forms caustics. 

Stretching his stiff legs, Persikov got up, returned to his laboratory, yawned, rubbed his permanently inflamed eyelids, sat 
down on the stool and looked into the microscope . . . With his right eye Persikov saw the cloudy white plate and blurred pale 
amoebas on it, but in the middle of the plate sat a coloured tendril, like a female curl. The coloured streak of light merely 
got in the way and indicated that the specimen was out of focus. The zoologist's long fingers had already tightened on the 
knob, when suddenly they trembled and let go ... He noticed that one particular ray in the coloured tendril stood out more 
vividly and boldly than the others . . . This strip of red was teeming with life. The old amoebas were forming pseudopodia in 
a desperate effort to reach the red strip, and when they did they came to life, as if by magic. They split into two in the ray, 
and each of the parts became a new, fresh organism in a couple of seconds. 

Mikhail Bulgakov, The Fateful Eggs, 1924. 

1 Introduction 

Caustics are the envelopes of the light rays. At the caustic points, intensity of light is singularly large, 
causing different physical phenomena (such as the one observed by Professor Persikov). Mathematically, 
caustics could be characterized as points where usual bounds on oscillatory integrals are no longer valid. In 
this paper we are going to consider how this concept applies to Fourier integral operators. This question 
becomes interesting in view of a recent paper |JMR00j on dissipative semilinear oscillations, where the L q 
estimates on oscillatory integrals with caustics played the central role. Our goal is to investigate how the 
regularity properties of Fourier integral operators are affected by the presence of caustics. We will show, in 
particular, that for q away from a certain neighborhood of q = p' the L p — > L q bounds on Fourier integral 
operators are caustic-insensitive. 

Oscillatory integrals with caustics have enjoyed much attention. The classical references are jAC ZV88 
and |Dui74| . The asymptotics of oscillatory integrals near caustics were derived in |Lud66| and |GS77| . 

Let us mention previously known estimates on Fourier integral operators. The L 2 estimates on Fourier 
integral operators were considered by Hormander H6r7lj. The LP — > L q and LP — + LP estimates on Fourier 
integral operators in the context of strictly hyperbolic equations with constant coefficients were addressed 
in |Str70| . |Lit73| . |Bre75| . and |Sug94| |Sug96l |Sug98| . L p -> L p ' estimates for certain hyperbolic equations 
with smooth coefficients and applications to the existence and uniqueness results for semi-linear hyperbolic 
equations are in ,Bro77i. The L p — > LP estimates were derived by Seeger, Sogge, and Stein J533S91 . For 
more information on regularity properties of generalized Radon transforms and Fourier integral operators 
associated to local graphs and to degenerate canonical relations see the reviews |GSW00| 1GS02 . 

•This work was supported in part by the NSF under Grants No. 0296036 and 0200880 
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We first recall some background about caustics of oscillatory integrals. Let us consider an oscillatory 
integral 

u T {x)=T k ' 2 [ e lT ^ x - a) a(x,T,a)d k a, a£R k , r > 0. (1.1) 

We assume that t/j is a smooth function and that a 6 S d is a symbol of order c? in r, compactly supported 
in a and x. If there are no critical points of the map a i— > V'C 3 '; a )i so that V'q 7^ everywhere in an open 
neighborhood of the support of a(x, t, a), then the repeated integration by parts shows that |u T (x)| = o(t~ n ), 
for any N > 0. If there are non-degenerate critical points, where ip' a = but detV'aa 0j then the method 
of stationary phase shows that |u r (x)| = 0(r d ). If we also assume that rank^" Q > n (when k > n), then one 
can readily show that ||u T (x)||i2 = 0(r d ). It follows that as long as all critical points are non-degenerate, 
u T (x) £ L q (R n ), 2 < q < oo, with the norms bounded uniformly in r £ (0, oo). 

If there are degenerate critical points, known as caustics, then ||u T (a;)||L° o is no longer uniformly bounded. 
The order of a caustic k is defined as the infimum of k' so that ||it T (a;)||x,oo — 0(t k ). For example, 
ip(x, a) — a 3 + xa corresponds to the fold (A 2 ), with k = 1/6; ip(x, a) = a 4 + X\0? + x-i corresponds to the 
cusp (A3), with k = 1/4. For more details, see |AGZV88] . |Dui74j . |Dui96| . At the same time, it was shown 
in |JMR00| that there exists q c > 2 such that the L q estimates for 2 < q < q c are still bounded uniformly in r. 
(This information was used to deduce that the singularities of solutions to dissipativc scmilincar equations, 
like Du = ii\u\ p ~ 1 , are absorbed at the caustic if p is larger than certain critical value; for generic caustics, 
one needs p > 3 for such an absorption to take place.) 

Now we turn to Fourier integral operators. Let X and Y be two smooth manifolds (without boundary). 
A Fourier integral operator $ : C^°(Y) — > @'{X) can be defined (locally) by 

$u(x)= f e i ^ x ' e ^a{x,e,y)u{y)dedy, (1.2) 

JR N xY 

where a is a symbol of order d and is a non-degenerate phase function. We write J € I^^X, Y, C), where 
the order of the operator is defined by \x = d + ^ — dlraX + dlmY an( j Q [ s the associated canonical relation. 
We will always assume that 

dim X = dim Y = n 

and that the symbol a is compactly supported in X x Y. 

Let us consider L 1 — * L°° estimates on J. From (|1.2J) one can see that 

d:L 1 (Y)^L°°(X) (1.3) 

if d + N < (which is equivalent with /1 < — (n + N)/2). The smaller the minimal number of oscillatory 
variables is, the better L 1 — * L°° regularity properties # possesses. As we know from |Hor71| . the minimal 
number of oscillatory variables is equal to iV m j n = 2n — r, where r is the minimal value of the rank of the 
projection ttxxy from C onto X xY: 

r = min rank ditx x Y ■ 

We define the caustic set of the canonical relation as a subset £(C) of C where the rank of dirxxY is not 
maximal: 

E(C) = {peC: rankd7rxxy| s(C) < 2n - 1}, 

so that outside of S(C) the number of oscillatory variables of a Fourier integral operator J £ I ti (X,Y, C) 
could be reduced to N = 1. Let J £ I^{X, Y, C), and let A > 1, be its Littlewood-Paley decomposition. 
Similarly to |Dui96| , we will say that k is the highest order of caustics of C if it is the infimum of numbers 
k' such that the Schwartz kernel of which is an oscillatory function of order /i, is bounded by 0(A' t+K ), 
uniformly in x and y. It follows that for the action 

5: L\Y)^L°°(X) 
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to be continuous we need to have 3 S 1^ with // < — ^yl — k. Thus, in the presence of caustics, the L 1 — > L°° 
estimates deteriorate. On the other hand, if we assume that C is a local graph, the mappings 

if M < 0, (1.4) 

n — 1 , 
if M< (1.5) 

if M<-f, (1.6) 

are continuous, independently of the presence of caustics. (The estimate (11.4(1 is the classical L? bound on 
Fourier integral operators, ((1.5(1 is proved in SSSjH], and ((1.6(1 follows from the f) 1 — * L°° estimate on #3* , 
which is a pseudodifferential operator of order 2/i.) The LP — > L 9 estimates for 1 < p < gr < 2 (obtained by 
interpolation of ((1.4(1 . 1(1.5(1 . and ((1.6(1 ^1 are also caustic-insensitive. By duality considerations, the same is 
true for 2 < p < q < oo. We are going to show that for l<p<2<<7<oo away from a certain neighborhood 
of the line q — p' the LP — » L q estimates are also caustic-insensitive. In this paper, we only consider the 
situation when rankditxxY > 2n — 2. 



3: 


L\Y) - 




3: 




,L\X) 


3: 




- L\X) 



Let us give a short account of our methods. Let C C T*(X)\0 x T*(Y)\0 be a canonical relation which is 
a local graph. Let 3 = ^ A 3a + 3o, A = 2 fe , fc G N be a Littlewood-Paley decomposition of 3 G I^iX, Y, C). 
If the caustic set £(C) is empty (rankd7TxxY = 2rt — 1 everywhere), then, representing 3 with fleR 1 , one 
easily checks that 

HA^^aIIl^l- <C, (1.7) 

uniformly in A. Now let S(C) 7^ 0, and assume that r&nkdirxxY = 2n — 2 at £(C). Choosing local 
coordinates a on the unit sphere in R N , we introduce a function 

V = det(\d\- 1 C ), l<i,j<N-l, (1.8) 

which measures the distance to the caustic set. This function is defined up to a nonzero factor, which depends 
on the local coordinates. We decompose 3a into E ct 3a, ct + 3a, nice, where a = 2~- J , j 6 N, with respect to 
the values of T>, so that the Schwartz kernel of 3a, a is localized to the set a/2 < \D\ < 2a near E(C), while 
the Schwartz kernel of 3a, nice is localized away from £(C). When approaching the caustic set, the L 1 — ► L°° 
estimates become worse: 

IIA^-^a^IIl^l- ~ct- 1/2 . (1.9) 

This is the optimal estimate which one expects from the application of the stationary phase method. On 
the other hand, the L 1 — > L 2 action of 3a, a improves near S(C): 

||A-"-*3 r A,<r|Ui^i»~^. (1.10) 

Similarly to the idea from |Tom79| . this estimate is essentially the "square root" of the estimate on the 
L 1 — > L°° action of A~ 2m ~"3a JS\,<t- The Schwartz kernel of A~ 2m ~™3^ .3a.ct is bounded uniformly in x, y, 
and A, and hence this operator is bounded from L 1 to L°° (uniformly in A). Moreover, the Schwartz kernel 
involves an inert integration in 6, and if T> vanishes of order m with respect to 9 (as for the caustics of the 
type A m+ i), then there is an improvement ||A -2 ' 1 ~"3a -3a,(t||l 1 ^l o0 ~ cr 1 /™ for small values of a, which 
leads to (fTTTOf) . 

The estimates (11.9(1 and 1(1. 1UH allow us to prove that there is some q c > 2 such that for 2 < q < q c the 
L 1 — > L q regularity of 3a is not affected by caustics. In essence, this situation is expressed by the following 
obvious lemma: 

Lemma 1.1 Let Bt, < t < 1, be a family of complete Banach spaces and that B t C Bf , \\ ■ \\b ± > || ■ ||s t , 
if t >t' . Moreover, assume that for any £ G B% C Bq, 

Let £j G B\, j G N, be a sequence such that ||£j||s < aP , < a < 1, and that \\£j\\B! < b J , b > 1. Then 
Ej-gnO converges in B t for <t <t c = _ hl ^ + ° lnb - 
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This already allows us to calculate the critical value q c . The estimate (II. 9f) (mapping to L°°) blows up as 
cr -1 / 2 , while the estimate (|1.1(J|1 (mapping to L?) improves as a 1 ^ 2 " 1 ^. Interpolation shows that the bound 
on the L 1 — » L q mapping of \~ ,J '~( n+1 M 2+1 / q $\ lCr behaves as o- - V^Ci+V" 1 )/^ which improves for small a if 
q < q c = 2 + 2/m, so that the mapping 

A -M-(n+l)/2+l/^ A . L l {Y) ^ L g (x); 2<q<q c , 

is bounded uniformly in A and is not affected by the caustics. These estimates could be interpolated with 
the (caustic- insensitive) L 2 — > L 2 estimates on For l<p<2<g<oo, Littlewood-Paley theory implies 
that the mapping $ : LP a — > L q is continuous if the mappings A"Sa : LP — > L 9 are bounded uniformly in A. 
This, together with duality considerations, yields the range of p and q such that the L p — > estimates are 
caustic- insensitive. 

Remark 1.2 This situation is similar to the LP — ► L p regularity of Fourier integral operators associated 
to degenerate canonical relations, when the projections from the canonical relation are allowed to have 
singularities. For example, as we know from |MT85j . if both projections C -> T*(X), C -> T*(Y) have at 
most Whitney fold singularities, then the operator # £ I^(X, Y, C) loses 1/6 of a derivative in the Sobolev 
spaces: $ : H a (Y) — > iJ" Q ~ M ~e(A), but according to |SS94j the L p — > L p regularity of such an operator 
for p £ (1,3/2) U (3, oo) is the same as for operators associated to local graphs |SSS91| : # : L^(Y) — > 
-^o!-^-(n-i)|(5 |(^0> = p — I" Similar estimates on operators with one-sided Whitney folds were derived 
in [UC03 . This time, one uses the Phong-Stein decomposition |PS91| of # with respect to the distance to 
the critical variety £ (where the projections from the canonical relation become singular). This distance is 
measured by the function 



h = |6>f " n det 



A" AJ> 



which is proportional to the determinants of the Jacobi matrices of projections from C. (The factor in the 
definition of h is chosen so that h is homogeneous of degree in 9.) We decompose $ = J2h^ h + Ssmooth, 
where h = 2~ J , j £ N. The operator is obtained from $ by localizing its integral kernel to the variety 
where ft/2 < \h\ < 2H, and the projections from C have no singularities on the support of the integral kernel 
of Ssmooth- The main observation is that while near E the Sobolev estimates become worse, HStiIIi, 2 — >£ a ~ 
hr 1 ! 2 , the Hardy space to L 1 estimates improve due to smaller size of the support of the integral kernel: 



Caustics of Lagrangian distributions are discussed in Section [2J The main results (Theorems 13.51 and 
13.11(1 are stated in Section [3J The proof of L p — > L q estimates is in Section QJ The sharp f) 1 — * L q estimates 
are proved in Section |SJ We apply our results to the estimates on the half- wave operator in Section [S] 

The consistency of definition (|1.8fl of the distance V to the caustic set is proved in Appendix [X] The 
technical lemma (f) 1 — > L°° bounds on pieces) which allows us to obtain t) 1 — > L q estimates is proved in 
Appendix [B] 

2 Caustics of Lagrangian distributions 
2.1 Symbols 

We will use the class of classical (polyhomogeneous) symbols, in the sense of [Hot94|. 



Definition 2.1 A smooth function a(x,9) on X x R N is called a symbol of order d in 9 if for any multi- 

?N 
J + 



indices a £ Z" and (3 £ Til 



\d%d$a(x,6)\ <G a ^(l + \9\) d ' m , for all (x, 9) £ X x R w , 
where \/3\ = j3\ + . . . + Pn- 
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We denote the class of symbols of order d by S d (X x R N ) or simply by S . 

The class of classical (or polyhomogeneous) symbols S^Xxt*) consists of symbols a(x,9) £ S d (X xM. N ) 
that satisfy an asymptotic development of the form 

oo 
3=0 

where aj are smooth functions on X x R w positively homogeneous of degree d — j for \0\ > 1; 

a 3 (x,T0) = T d - j aj(x,6) if \6\>1, r > 1. (2.2) 
The asymptotic development (I2.1|l means that we have 

fc-i 

a(x,6)-J2 a i( x > d ) =0(\6\ d - k ) for |0| > 1 (2.3) 

3=0 

and similar estimates for the derivatives. 
2.2 Oscillatory functions 

Let X be a C°° manifold and A C T*(X) a C°° Lagrangian submanifold. We say that ip £ C°°(A x M fe ) 
parametrizes A (locally) if 

d( x ,a)da'4' nas rank /c when d a ip = (2-4) 

and A is locally given by 

A^ = {(x, dxi)(x, a)) : d a ip(x, a) = 0}. (2.5) 

Definition 2.2 Let A be a C°° Lagrange submanifold in T*(X). An oscillatory function u(x,t) of order fi 
defined by A is a locally finite (in X) sum of integrals of the form 

I(x,t)= [ e tT ^ x ' T ' a) b(x,T,a)da, 



____ W~f" 

where a = {a\ , . . . , otk), k G N, %jj satisfies \2.4\) , A^ is a piece of A and b(x, r, a) G S cl 2 (a classical symbol 
of order \i + -| in t) vanishes for a outside a fixed compact set in M fe . 

Definition 2.3 (Duistermaat [Dui96j) Let i : A — > T*[X) be an immersed Lagrange manifold in T*(X). 
The caustic c(A) of A is the projection into X of the set S(A) of points in i{A) where i is not transversal 
to the fibers. At each point x £ X the order of the caustic is defined as the infimum k(xq) of the numbers 
K 1 such that u(x,t) — 0(r A ' +K ) for r — > oo ; uniformly for x in a neighborhood of x , for any oscillatory 
function u of order \x defined by A. 

Of course k(xq) = for xo € 7r(A)\c(A) and k(xq) < k/2 where k is the maximum of the dimensions of the 
intersections TV X0) £)(A) fl Tr XOt g\ (fiber) where (xo,^) S A. 

2.3 Lagrangian distributions 

We also need to define caustics of the conic Lagrangian submanifolds. Let A be a C°° manifold and 
A C T*{X) be a conic C°° Lagrangian submanifold. We say that <j> £ C°°(A x R N ) parametrizes A (locally) 
if 

d( x ,9)de4> has rank N when d$<fi = (2-6) 

and A is locally given by 

A = {(x, d x( j>{x, 9)) : d e <f>(x, 6) = 0}. (2.7) 
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Definition 2.4 (Hormander) Let A be a C°° conic Lagrangian submanifold in T*(X). A Lagrangian 
distribution u(x) of order p defined by A is a locally finite (in X) sum of integrals of the form 

where 9 G WL N , <f> satisfies $2.b)) . A^ is a piece of A and 

a(x,6) e S^^+^iX x R N ). 

We pick a smooth function p e Cg°([-2,2]), p > 0, a] = 1. Define (3{t) = p{t) - p{2t) for t > 0, 
(3 = for i < 0, so that /3 € Cq°([|,2]). We introduce the Littlewood-Paley decomposition of u(x): 

u x (x)= J e^ {xfi) [3{\9\/\)a(x,9)d9. (2.8) 

Definition 2.5 Lei A be a conic Lagrangian manifold in T*{X). The caustic c(A) of A is the projection 
into X of the set £(A) of points in A where rankdTTx < dim A" — 1. At each point xo € X the order of the 
caustic is defined as the infimum k(xq) of the numbers k' such that u\(x) = 0(A Al+d ' 1 4 X +k+ K j for A — > oo, 
uniformly for x in a neighborhood of Xq, for any Lagrangian distribution u of order fi defined by A. 

Definition 2.6 We say that A has a caustic of corank k > 1 at a point po £ S(A) if 

rank(i7rv| = dim A" — 1 — k. 

Lemma 2.7 Let A be a smooth closed conic Lagrangian submanifold ofT*(X)\0. Let <j>(x, 9) e C°°(XxR N ) 
be a smooth non-degenerate phase function which parametrizes A: 

A = {(x,d x cb(x,9)) : do<Kx,6) = 0}. 

Let a = {oti}, 1 < i < N — 1, be local coordinates on the unit sphere E> N_1 . We use (A, a) € M+ x § w_1 as 
local coordinates in R N . Then T> = detij(A -1 0^. Q . | A ), 1 < i, j < N — 1, is a smooth function on A defined 
up to a nonzero factor: 

Der(A)/C™(A). 

This statement is intuitively clear, since T> vanishes precisely on the caustic set E(A) where the rank of the 
projection from A onto X is smaller than dim X — 1. Still, since we need to know that the order of vanishing 
of T> at S(A) in particular directions does not depend on the number of oscillatory variables and the choice 
of local coordinates, we will give a detailed argument in Appendix 1X1 

Definition 2.8 We say that the caustic at a point po € S(A) is simple if it is of corank k — 1, so that 

rankrf7ry| = dim AT — 1 — k = dim AT — 2, 

and if d( x j ) V\ m ± 0. 

Definition 2.9 We say that the simple caustic at a point po £ S(A) is of index m 6 N if m is the smallest 
integer so that there exists a vector field V € C°°(r(T(A))), € hevdirx, such that 

V m V( Po ) ? 0, 

where 

x> = det(|0rV' ), i < *, j < jv - 1. 

Remark 2.10 Consistency of this definition (independence of the choice of the phase function <f> that 
parametrizes A and independence of the choice of local coordinates a on S 1 ^^ 1 ) follows from Lemma \2.T\ 
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Example 2.11 Let 9 G R 2 and A = Then 6/\9\ £ S 1 . Denote a local coordinate on S 1 6ja. Consider 
the phase function <p(x, 9) = \6\<$>(x, a), with 

®(x, a) = a m+2 + cciQ! m + . . . + x m a + x m+ \. 

This is the model example of a caustic of the type A m+ i. The corresponding Lagrangian is A — {x, d<p(x, 9) : 4>'e = 
0}, which can be written as 

A = {x, \d x $(x, a) : $(x, a) = 0, & a (x, a) = 0}. 

The Lagrangian could be parametrized by (x',X,a), where x' — (xx, . . . ,x m —i). At A, one can express x rn 
and x m+1 as x m = X m (x',a), x m+1 = X m+1 (x', a). 

The kernel of the map it : A — > X always contains the tangent vector d/d\. On the caustic set 

c(A) = {(x',a) G A : V(x',a) = ^„ a (x', x m (x' , a), x m+1 (x' , a), a) = 0} 

one also has dj da £ ker<i7r. Since 1 < dimkerd7T < 2 and dT> — d(^>'^ a \ A ) ^ 0, the caustics are simple in 
the sense of Definition \2.fA 

Consider the vector field V = d/da G C°°(T(TA)), V\ c(A) G kerdTr. Since 

V m V (x', a) = 8% (K a I a ) = W , *m (x' , a) , x m+ x (x' , a) , a) + 0, 

one concludes that the caustic is of index (at most) m. 

Remark 2.12 While caustics of the type A m+ i, m > 1, correspond to simple caustics of index ra, the 
converse is not necessarily true, except when m = 1 and 2. 

Let us show how to prove that simple caustics of index m = 1 and 2 necessarily correspond to caustics 
of the type Ai and A3, respectively. We first reduce the number of oscillatory variables to N = 2, denote by 
a the local coordinate on S 1 , and define 

*(,x,a) = <l>(x,6/\6\)=<l>(x,e)/\9\. 

It suffices to notice that </> has the caustic of the type A m+1 at the point (xq, ao) if 9^$(xo, ao) = 0, j < m+2, 
d™ +2 §(x, a) ^ 0, and if the differentials 

d$(x, a), d$' Q (a;,a), d$^ m) (x, a) 

are linearly independent. In the case m = 1, the linear independence of d$ and d$^, follows from the 
non-degeneracy assumption on <f) (the differentials dcj>g- are linearly independent). 

To settle the case m = 2, we additionally need to check that d^'^ a is linearly independent of d$ and 
d$' a . We only need to notice that two latter differentials vanish identically on TA, while the differential 
rf$JJ„ = dT> was assumed to be different from zero (see Definition 12. 8|l . 

3 Main results 

Consider a Fourier integral operator 

du(x) = f e^ x ' e ^a{x,9,y)u(y)d9dy, (3.1) 

JR N xY 

where X and Y are two smooth manifolds (unless stated otherwise, we assume that dimX = dim!" = n), 
a(x, 9, y) € Sf y (X x x Y) is a symbol of order d in 0, compactly supported in X x F. (We restrict the 
consideration to classical (polyhomogeneous) symbols, denoted by S c \.) The function is a non-degenerate 
phase: the differentials d(a;,e,y)^e-) 1 — J — ^> are linearly independent, so that 

S^={(x,9,y):cf>' e (x,9,y)=0} 
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is a smooth submanifold of X x M. N x Y of dimension dimJT + dim Y". We write # € I^(X, Y, C), where the 
order of the operator is defined by /i — d + ^ — dlmX + dlmY ana ; q [ s the associated canonical relation: 

C = {(x, 0, y)), (y, -d v <j)(x, 6, y)) : cp' (x, 9, y) = 0}. 

According to |Hor71| , the minimal number of oscillatory variables is equal to 

TV m in = dim X + dim Y — min rank dux x y ■ 
If C has non-empty caustic set S(C), then iV m i n > 1. 

Definition 3.1 We will say that the canonical relation C C T*(X)\0 x T*(Y)\0 has a caustic at a point 
po = ((xo, £o)> (j/0) "Ho)) G C if the twisted canonical relation C = {((x, £), (y, — 77)) : ((x,£),(y,r))) £ C}, 
which is a conic Lagrangian submanifold of 'T*{X x V)\0, has a caustic at a point ((xq, £o), (yo, — ?7o)) C C. 
We will not distinguish the caustics of C and C. 

The following result is an immediate consequence of Definitions 12.51 and 13. II 

Theorem 3.2 Let X , Y be two smooth manifolds (possibly of different dimension), and let C C T*(X)\0 x 
T*(Y")\0 be a smooth canonical relation. Let the Fourier integral operator $ G L tJ '(X, Y, C) have its symbol 
compactly supported in X x Y . (The symbol of ^ does not have to be poly homogeneous.) If C has caustics 
of order at most k, then 

i , ,* / \ dim X + dim Y 1 , 
5: L 1 (y)^i°°(X) if ^< 1 --«. (3.2) 



Further, assume that dim X = dim Y~ = n and i/iai C is a local graph. Then 

1 

P 2 



: Ll +{n+1+2K)Sp (Y) -> 1< p < 2, ^ = 1 - 1. (3.3) 



Proof. The first part of the theorem follows from the trivial estimate 

HffAulU- < CA^ + dimX t dimy + H«|| M || Ll . (3.4) 
For the second part, we interpolate l|3.4|) with ||5^a || i 2 — »i 2 < CA P and apply Littlewood-Paley theory. □ 

Definition 3.3 For our convenience, we introduce the map 

(•,-) f : (P,q)^(p,q) i = (l/pA/q)- (3.5) 

Definition 3.4 We define 

2 2 

Pm = 2 — , g m = p' m = 2 H 

m + 2 to 

and introduce the following regions in (1/p, 1/q) -plane (see FiaureU^: 

2l m is the open triangle with the vertices at (2,2y, (1, ly , and (l,q m y. 
S m is the open triangle with the vertices at (oo,oo)t. (2,2)^, and (p m ,ocy. 
(£ m is tte open convex hull of points (1,oo)t, (l t Qm) > (2,2)t, and (p m ,cxD)^. 

The following is our main result: 

Theorem 3.5 Let X, Y be two smooth manifolds, dimX = dimY" = n, and let C C T*(X)\0 x T*(Y)\0 
be a smooth canonical relation which is a local graph. Assume that C has only simple caustics of index at 
most m, m € N. Let 3" £ /^(X, Y", C) /lave t/ie classical (polyhomogeneous) symbol with compact support in 
X x Y. Then for (p, g)t G 2l m U Q3 m i/ie i p -> i« estimates on $ are caustics-insensitive. 
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p = oo 



Pm = 2 



m + 2 



Qm = P' m = 2 + — 

m 



I 

p = 2 P = Pm p= 1 p 

Figure 1: Regions 2t m <B m , and £ m in (l/p, l/g)-plane. 



Precisely, 



whe 



3-- Ll +nSp+Sq (Y) ^ L«(X), (p^jtgSU, 



<5d = 



1 1 



p 2 '29 
For (p, 9)^ G Cm, fie estimates depend on the order of the caustic, given by K = \ 



1 . 

m+2 ' 



$ : L M+^p + (^+^)(l/2+K) + (5 9 -5 P )( r ) G Cm, 9 < P, 

Sq){1/2+K)+{Sp _ 5q) (Y) ^ L«(X), (p,q^eC m , q>p'. 



^ ' ^ ti+nS q + (S p +u q 



(3.6) 
(3.7) 



(3.8) 
(3.9) 



Remark 3.6 We restrict the consideration to the class of classical symbols C Sf Q in order to simplify 
the proof of Lemma \4-'J\ 



Remark 3.7 Note that k 



\ ~ m+2 * s the order of a caustic of the type A m+ \. 



Remark 3.8 We need the assumption that C is a local graph in order to interpolate with the L 2 -based 
Sobolev estimates on 

ff: Ll(Y)^L 2 (X). (3.10) 

The argument could immediately be adapted to the case when the projections from C have singularities, as 
long as C — > X and C — > Y are assumed to be submersions. In this case, one only needs to modify QZS.lUp . 
taking into account the loss of derivatives due to singularities of the projections; see \GS02L 



Remark 3.9 The sharp estimates on the line p = q follow from )SSS911f : 

5: t)l + ^i(Y)^L\X). 

This map can be interpolated with the continuous L 2 ^ — > L 2 action. ( Generalization for operators with 
degenerate canonical relations is obtained in \CC03f . ) 

Remark 3.10 On the line segments ((2, 2)\ (1, q m ) )> ((2, 2y , (p m , oo)^) 7 and on the lines p = 1 and q = 00 
the stated estimates hold with the loss of e > 0. 
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In particular cases, we also have sharp t) 1 — > L q and LP — > BMO estimates, as stated in the next theorem. 

Theorem 3.11 Let X, Y be two smooth manifolds, dimX = dim Y = n, and let C C T*(X)\0 x T*(Y)\0 
&e a smooth canonical relation such that both C — > X and C — > Y are submersions. Assume that C /las oniy 
caustics of the type A m+ \ with m = 1 or 2. Let $ S L^(X, Y. C) ftawe i/ie classical (polyhomogeneous) symbol 
compactly supported in X x Y . Then 

5 : ^ . » , x , ( y ) «™ < g < oo. (3.12) 

Remark 3.12 In this theorem, we do not need C to be a local graph. 

Remark 3.13 L' P (Y) — > BMO(X) estimates on $ for 1 < p < 2 are obtained by duality. Other LP — ► L 9 
estimates can be obtained by interpolation with L 2 Sobolev estimates. 

4 Microlocal techniques: decompositions and interpolations 

In this section, we prove Theorem 13. 51 



4.1 Dyadic decompositions 

We pick a smooth function p € C* °°([-2, 2]), p > 0, a] = 1. Define /3 € Cq°([|, 2]) by /3(t) = p(t) - p(2t) 
for i > 0, (3 = for t < 0. The functions p and /3 define dyadic partition of unity: for any tel. 

x;e^ ±2 " j */2)+p(i*d=i. 

± jeN 

We use the partition of unity which is the Littlewood-Paley decomposition with respect to the magnitude of 
\9\ and the dyadic decomposition with respect to the distance V from S(C): 

1=1 ]T fi{2- l \e\)+p{\e\)\ [^^/3(±2^)+p(2^|P|) + (l-p(2|P|)) 

\A=2 ! ,ieN / \ ± j=l 

We define 

$x,±*u(x)= f e^ (x < e ^P(±V(x,8,y)/o)[3(\9\/\)a(x,6,y)u(y)dddy, (4.1) 

JR N xY 



We also define 



d\,*u(x)= / e^ x ' e ^p(V(x,0,y)/a)/3(\e\/\)a(x i e,y)u(y)d0dy. (4.2) 

Jl«xY 

ff S mooth«(a;) = / e^^V(|0|MM,yMy)d^?/ (4.3) 
and 

$niceu(x) = f e***^ [(1 - p(\6\)) (1 - p(2\V\)] a(x, 9, y)u{y) dOdy. (4.4) 



'xY 



There is a decomposition 



5 = X!X! X! ^ A .±°- + X^^oM + Smooth + Snicc, (4.5) 
± A cr>cr (A) A 

where both A and o run over powers of 2: 

A = 2 l , ZeN, o = 2-l, i< j<j (A) = [[log 2 A^N]]. 
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We set 00(A) = 2~^ x \ so that 

a (X) w A"^. (4.6) 

We use the symbol "w" to indicate that the quantities differ at most by a factor of 2. 

The operator ^smooth is infinitely smoothing and can be discarded. Since there are no caustics on the 
support of 1 — p(2|2?|), the operator 3m CO can &l so be discarded. The estimates on operators $\.± a are the 
same independent of the sign, and the treatment is the same; we will only consider the "+"-case. 

4.2 L 1 — > L°° estimates 

Proposition 4.1 Let C and $ S i 7 * (X, Y, C) be as in Theorem I.V.,51 and let $\ tCr , 3A,<r be given by 1)4. 
(TO) . T/ien 

II^Hl^l- <CA^+^(t-5, (4.7) 

IIMn-i- + II5a, ct ||l^l- < CA^+^a^. (4.8) 

Remark 4.2 TTie uahte (To (A) w A~ m + 2 m j^.6| ) is chosen so that the estimates \4- j4-8{ ) coincide at 
a = (To(A). 

Proof. We use the representation of 5 with the minimal possible number of oscillatory variables, N = 2. 
Then $\ a could be written as 



e 

xS JY 



j+ { *< T ' a,v) pWL(*, 1, «. v)/<r)P(r/\)o(x, r, a, p)u(y) rdr da dy, (4.9) 
where a(x, t, a, y) is a classical symbol of order d = fj, + § — 1. For H4.70 . we need the bound 
e ^(^.«.»)p(^ a ( X) 1, «, y)/a)P(T/\)a(x, r, a, y) rdr da 



xS 



< CX^+^a^i, (4.10) 



uniformly in x and y from a small open neighborhood in X x Y and for all A > 1 and er < 1, a > A - ™+ 2 . 
For simplicity we assume that /i = — n/2, so that ra(x, r, a, y) is a symbol of order zero, which we denote 
by 6(x, r, a, y). This classical symbol has the development 

b(x,r,a,y) ~ b Q (x,a,y) + ^2b j (x,a,y)T~ : > , t > 1. (4.11) 

Denote 

= A~ 1/2 /" dT /" ^ e **(*,r,a,v) 6(ajT)as ^ (T/AW ^ (a;jl)a}y)/£r)) (4 . 12) 
where C § denotes a-support of b(x, t, a, y). Substituting t = Az, we rewrite I\. a {x, y) as 

I X)(T (x,y) = X 1 / 2 [ dz f dae iXz ^ a ^b(x,Xz,a,y)f3(z)P(^ a (x,l,a,y)/a). (4.13) 
For H4.10f> . we need the bound 

\hA^y)\ <c<j- xi \ (4.i4) 

valid for all A > 1, \- m /( m + 2 ) < a <1, and with C independent of A and a. 
Lemma 4.3 

\h,a(x,y)\ <C<J- 1/2 , (4.15) 
with C < oo independent on X > 1 and < er < 1. 
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We need this estimate to be uniform in A and a simultaneously. Similar estimates were considered in 
CdV77 and in many other papers. The result is known to be optimal, but is not proved in the whole 
generality in higher dimensions. For the sake of completeness, we give our own proof for the case we are 
interested in. 

Proof. There is a trivial bound \I\ y(J {x, y)\ < A 1 / 2 , due to the compact support of the integrand. This settles 
the case < a < A ; from now on, we assume that A -1 < a < 1. 
Denote b'(x, t, a, y) — b(x, t, a, y) — bo(x, a, y) 6 ^ci 1 , and let 

I' x , a {x,y) = A" 1 / 2 f dT f dae lXz ^ 1 -^b'(x,Xz,a,y)f3(z)P(^ a (x,l,a,y)/a). (4.16) 
Jm. Jk 

Since \b(x,Xz,a,y)[3(z)\ < CA _1 , uniformly in x, a, A > 1, 1/2 < z < 2, and y, there is an easy bound 

\l' x<(7 (x,y)\ < A 1 / 2 f f ^{xAz^a^^dzdaKCX- 1 ' 2 <C. 
Ji/2 Jk 

Thus, we only need to consider the bound on l|4.13[) with bo(x, &i y) instead of b(x, Az, a, y): 

Il a (x,y) = \ 1 ' 2 [ [ e iX ^ 1 ^h (x,Xz,a > y)[3(z)f3(^ a (x,l,a > y)/a)dzda. (4.17) 
Jn Jk 

Denoting by (3 the Fourier transform of /3(z), we rewrite I x a as 

Il„(x, y) = A 1 / 2 f b (x, a, y)0(X<j>(x, 1, a, y))p(K a (», 1, a, y)M da. (4.18) 
Jk 

The statement of the lemma follows from the bound 

\li a {x,y)\<Co-V\ (4.19) 

which is uniform in A, A -1 < a < 1, x, and y. We will prove this bound in the next lemma. 

This completes the proof of Lemma f4. 31 □ 



Lemma 4.4 Assume that 0Q Q (a) vanishes at most of order m on a compact set K C 
/// e L 1 (K)n J L 00 (]R) and [3 e C£°([l/2, 2]), then 

A 1 / 2 / fiXmmKai^/^da^ilfh. + lfh^Ca- 1 / 2 , 
Jk 

uniformly in X > 1 and < S < 1. 

Essentially, we are proving the following sublcvel set estimate: 



(4.20) 



{aeK: |0(a) - 7 |<A-\ > a 1 / 2 } 



< 



C 



Xa 



(uniformly in 7 € K). 



1. We 



Proof. The proof of this estimate is simple, so we can give it in detail. Let p € C^°([— 2, 2]), p\ [ _ 1 „ 
use the partition 

1 = p(ff/ a VA^) + (l - pfay/Xfi)) 
to rewrite H4.20JI as a sum of two terms, 

^ [ fM)pW a V*M0(4>'L.M da + xi f f(Xcf>) (l - ptf a y/X/^)) 0(^Ja) da, (4.21) 



K 



K 



which we analyze separately. 
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The first term in (|4.21|) is bounded by 

A%|U~ / p{<t>' a {a)V^WLM)l<r)da 



K 



<A'||/|| 



G 



^infl^LI 



< CV 



l/ll 



since inf |</>^ Q | > cr/2 on the support of the integrand. The value of C depends on the bound on the number 
of roots of 4>' a (a) = c (this number is bounded uniformly in c due to the finite type assumption: <j)'^ a vanishes 
of order at most m). 

The second term is bounded by 

C 



A 1 / 2 I /(A#a))(l - p^MVW)) da < X^f\\ L1J —^ < Co-^\\f\\ L1 



K 



since inf \(j)' a \ > \J a jX on the support of the integrand. Again, we need to mention that the number of roots 
of 4>a a ( a ) = c is bounded uniformly in c due to the finite type assumption. 

This proves Lemma R~4l □ 



Remark 4.5 The maximal order of vanishing, m £ N, does not appear in the above lemma. The statement 
of the lemma is also true without the finite type assumption if we require that <f> is real analytic, or, more 
generally, if we require that <f> S C°°(R) and that (f>" is "finitely oscillating" on K C K: 

Number of connected components of the set {a £ K : (f>"(a) — c} is bounded uniformly in c G R. 

This assumption holds for any real analytic function, but does not hold for all smooth functions; an example 
of a smooth function which is "infinitely oscillating" on [—1, 1] is er x l x sin(l/x). 

This finishes the proof of (|4.7(l . 

The proof of (|4.8() is similar but much more straightforward. One needs to use the following well-known 
lemma (see, e.g., |CCW99| 'l: 

Lemma 4.6 If f(a) vanishes at most of order m on [—2,2] and (3 G C^°([— 2,2]), then J R (3(f(a)/a)da is 
bounded by Co 1 !™ . 

This finishes the proof of Proposition ^. II □ 



4.3 L 1 — > L 2 estimates 

Lemma 4.7 Ul^xAW^ + Ul^xAW^ < CX 2 »+ n o^. 

Proof. Since we assume that C is a local graph (or at least that C — > Y is a submersion, as in Theorem 
13.11(1 . we can choose the phase function of the form 4>(x, 9,y) — x ■ 9 — S(6, y), with 9 £ M. N , N — n, where 
S(9,y) is homogeneous in 9 of degree 1. Then 9 and y can be used as the local coordinates on C. We can 
rewrite J in the form 

$ K *u(x) = I e^- e ' s( - e ^a(9, y)P{4>{9, y)/a)u(y) d9 dy. 

JM N xY 

3*3 £ I~ 2 ^ (X, X, A) is a Fourier integral operator with the phase S(9, z) — S(9,y), associated to the diagonal 
A c T*(Y) x T*(Y), and with N = n. 
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The Fourier integral operator 



e iW0,z)-s(0,y)) a ^ yHe , z)P(V{0, y)/a)0(V(0, z)/a)f3 2 (\9\/\)u(y) d9 dy 

R N xY 

has N = n oscillatory variables. (The number of oscillatory variables cannot be reduced since the rank of 
the matrix dg i dg j [S(9, z) — S(8, y)} is zero at y — z.) The order of its symbol is 2fi. This yields the bound 
const A 2ai+, V™ on the L 1 — ► L°° action, with the factor erm due to Lemma f4. 61 □ 

Remark 4.8 If C is a local graph, so that detij dg i d Vj S{9,y) ^ 0, then is a pseudodifferential operator. 
This lemma yields the following estimate: 

Proposition 4.9 Let C and $ 6 I^(X, Y, C) be as in Theorem \S.5l and 5a, o-, 3a, o- be given by (|4.1(l . (|4.2|l . 

Then 

UxAl^l* + WZxAl^ls < CX^a^. (4.22) 
4.4 IP — > L q estimates for 1 < p < 2 < q < oo 

Even when we can not prove the sharp fj 1 — > L q estimates (without the loss of e > 0), we still can prove the 
sharp LP — > L q estimates, for certain values of p and q which also satisfy l<p<2, 2<g<oo. The main 
tool is the Littlewood-Paley theory. 

We group the pieces Sa, a and 5a, a defined by (|4. 1|> . 14.2(1 into A-clusters: 

3\ = F\,<r+$\,<T W ( 4 - 23 ) 

ct=2-3, jeN, a>a (X) 

Let us consider this series in the norm of operators from Lj i+ „ +s (Y) to L q (X). 

Proposition 4.10 Let C and # € I^(X, Y, C) be as in Theorem \3.5l and let Sa be given by (|4.23() . Then 

Bxh^L, <C\^ +5 \ 2<q<q m , (4.24) 

WSxh^u < CA M+ * + *« +K ^^ , <? m < g < oo. (4.25) 

Proof. Interpolating L 1 — > i°° estimates from Proposition 14. II with L 1 — > L 2 estimates from Proposition 
we obtain: 



Lemma 4.11 



$A,ct||L1^L9 



< C\^ + ^ +s "a^- { ^ +1)s ". (4.26) 



Therefore, 1(4.23(1 is dominated by the geometric series. 

^ 2m ^ (m ^~ 1)^9 (equivalent with g < g m ) 7 then the geometric series is convergent, and hence is bounded 
uniformly in A. This proves ((4.24|) . 

If < (— + l)<5 g (equivalent with q > q m ), the series ((4.2311 is dominated by finitely many terms of the 
divergent geometric series: 

\\$xh^L« < C Yl X^+^o--^ +1 ^-^ (4.27) 

<r=2-J, jeN, <T><Tf,(A) 
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Taking into account that erg (A) = A m + 2 and that 



( ( JL + A $ - J_ ] = — ( 2 (i + m )A-i) (4.28) 
m + 2\\m J q 2m) 2(m + 2) V V ^ 9 ' y ' 



^(S g /8 qm " 1) = «^ 



<5« - 5, 



</.. 



we can rewrite (I4.27|) in a more convenient form: 

UxWl^l, < C\^ +5 « +K ^^ . (4.29) 
This proves (|4~TUI) . □ 

The estimates stated in Proposition 14. 1UI can be interpolated with the L 2 — > i 2 -estimates. If we assume 
that C is a local graph, then $\ ; ^ L 2 , and we obtain, for 1 < p < 2: 

\\$\\\l>^l* < CX» +n5p+s \ 2<q< 



1 - 4<5 p <5 9m ' 



h\\Lp^Li < C\ tJ - +nSp+( - 5p+Sq) ( 1/2+K ' >+{Sq ~ Sp \ <q <p'. 



According to Littlewood-Paley theory ( See93 , Lemma 2.1), 5 = J2\=2 l (sn^a nas the same LP — > 
regularity properties as long asl<p<2<<7< oo. This, together with the duality arguments, proves 
Theorem 13. 51 

5 Microlocal techniques: f) 1 — > L q estimates 

We are going to prove Theorem l3.11l which gives the substitute of the LP — > L q estimates for p = 1 (f) 1 — > 
estimates) and for g = oo (L p — > BMO estimates). 

5.1 f} 1 — > L°° estimates 

The following is the analogue of Proposition ^. II 

Proposition 5.1 Let C C T*(X)\0 x T*(Y)\0 be a smooth canonical relation such that C -> 1 is a 
submersion. Assume that C has only caustics of the type A m+ i with m = 1 or 2. Let $ G L^{X, Y, C) ftove 
£/ie polyhomogeneous symbol with compact support in X , Y , and let 3a, a- ? 3a, <t be given by l|4.1f> . I|4.2|l . Then, 
for any atom oq supported in the cube Q with side r, we have 

\\8\*a Q \\ L *> < CX^+^a-^ min(Ar, (Ar)" 1 ), (5.1) 

llffA.aOQik- + \\$\,*a Q \\ L ~ < CX^+^a^ min(Ar, (Ar)" 1 ). (5.2) 

Proof. The proof is similar to SS S91| . |CC03| . For the reader's convenience, we reproduce this proof in 
Appendix [BJ We require that a > A -1 / 2 (equivalent to m < 2) so that the localizations would not be too 
fine and the integration by parts from |SSS91j could be used verbatim. □ 

We group the pieces 3a, a into A-clusters as in (|4.23j) : 

3a = ^ 3a, o- + 5a,<t (a)- 

(7=2-3 , jgN, cr>cr (A) 

The estimates (|5.1|) and (|5.2() yield the following bounds on ||3AaQ||L°°: 
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Corollary 5.2 Assume that C — > X is a submersion and that C has only caustics of the type A m+ i with 
m = 1 or 2. Then 

llffAflQlU- <<7A^+^+"min(Ar, (Ar)- 1 ), K =i-_L_. (5.3) 
This allows us to conclude that 

5.2 f} 1 — > L 2 estimates 

Proposition l5.ll gives the sharp version of Proposition 14. II Now we are going to prove the sharp version of 
Proposition 14.91 

Lemma 5.3 Assume that C — > Y is a submersion and that C has only caustics of the type A m+ i with m = 1 
or 2. Then, for any atom oq supported in the cube Q with side r, we have 

1131,^,^(2 IU~ + \\Tx,Jx,*a Q \\ L ~ < C\ 2 » +n <r^ min(Ar, (Ar)- 1 ). (5.5) 

Proof. The proof is similar to the proof of Proposition 14.91 For the (j 1 — > L°° estimates, we can apply the 
usual machinery as long as miner w A - ™+ 2 is not smaller than A" 1 / 2 , that is, as long as m < 2. □ 

This lemma proves the following sharp version of Proposition 14.91 

Proposition 5.4 Let C C T*(X)\0 x T*(Y)\0 be a smooth canonical relation such that C — > Y" is a 
submersion. Assume that C /ias on/y caustics of the type A m +i with m = 1 or 2. Let ^ £ I^{X, Y, C) ftowe 
t/ie polyhomogeneous symbol with compact support in X , Y . and let 3a, a-, 3a, a be given by 14. ID . 14. 2f) . Then, 
/or any atom oq supported in the cube Q with side r, we have: 

Ux,„a Q \\ L 2 + U^a Q \\ L 2 < CA^+*a^min((Ar) 1 / 2 ,(Ar)- 1/2 ). (5.6) 

5.3 f} 1 — > L q estimates for small 2 < q < q m : cr-interpolation 

We group the pieces 3a, a- and 3a. a into c-clusters: 

3.- (5.7) 

A: cr>2cro(A) A: cr (A)<<7-<2cr (A) 

A=2 ! ,2eN A=2 ! ,2eN 

Then we have 

3 = ^ ' 3(T "H 3nice • 

Proposition 15 . II proves the following bound: 

Lemma 5.5 Assume that C X is a submersion and that C has only caustics of the type A m +i with 
m = 1 or 2 . Let 3^ &e given by (15.7(1 . TTien 

ll&rllo* ^Co-- 1 / 2 . (5.8) 

Proposition l5. 41 proves the following: 

Lemma 5.6 Assume that C — > V is a submersion and that C has only caustics of the type A m+ \ with m = 1 
or 2. Let 3a- oe given by H5.7[) . TTien 

Wfti ^L»<Ca^. (5.9) 
V+f 
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Corollary 5.7 Assume that both C — > A and C — > y are submersions and that C /las onZi/ caustics of the 
type A rn+ i with m — 1 or 2. Let 3v 6e given by J5.7| ). TTien i/ie interpolation of (|5.8|l and ()5.9|) gives 

ll^llfti „ < C«7^-(™ +1 ^, 2<g<oo. (5.10) 

The summation y^ (rT _ 9 _ 7 j^So- converges in f}^ + ri +5 — ► L q operator norm (where q > 2) if 



2m \m y 

which is equivalent to 2 < a < q m , q m = 2 + In this case, we conclude that 

5 : ^3+5, - i 9 , 2 < 9 < q m . (5.11) 
Note that the estimates (|5.11|l do not depend on the order of caustics. 

5.4 f} 1 — > L 9 estimates for g > g m : ^-interpolation 

In the case m < 2, we can derive the sharp f) 1 — > L ? estimates for q > q m . According to Proposition l5.il if 
m < 2 and if oq is an atom supported in the cube Q with side r, then 

HffA.aOQllic. + H^OQlli- < CA^ 2 * 2 ^ min(Ar, (Ar)- 1 ), (5.12) 

ll^agllL- < CA^+^+^A-V-' min(Ar, (Ar)" 1 ). (5.13) 

According to Corollarv l5.4l 

\\3x,*a Q \\ L 2 + \\Z x , a a Q \\ L 2 <CA^+t-^A^a^min((Ar) 1 / 2 ,(Ar)- 1 / 2 ). (5.14) 

We introduce a new parameter, u, for the values of A _K cr _ 5 (these values are bounded by 1 since a > 
cr (A) w A~ 2fi ). Let us group the operators Sa,o- into w-clusters $ u , u> — 2~ fc , k £ N, so that 

1? = Snice + ^ (5.15) 

L0=2- k , keti 

where 

a)<A" K (T" 1/2 <2a) w<A"' i CT" 1/2 <2w 
(T>2(7o(A) CT (A)<(T<2cro(A) 

K = 5 - A = 2', Z € N and tr = 2~ J ', j € N. 

Lemma 5.8 Assume that both C — > A and C — > Y are submersions and that C Zias onZi/ caustics of the 
type A m+ i with m = 1 or 2. Lei 3w be given by (|5.16|) . Then, for any atom ag supported in the cube Q wii/i 
side r, we ftave 

+1 < Clu, (5.17) 

HSwII^ „ B ^ < Cu- l ' m . (5.18) 

Corollary 5.9 Assume that both C — > A and C — > Y are submersions and that C has only caustics of the 
type A m+ i with m = 1 or 2. Let ^ u be given by (|5.16|l . Then 

WUW 2m+3 ^L«<Cu-^ +2 ^ + V 5 °. (5.19) 

m + 2 

The series X^=2-j jgn^ (considered in t)^ +2L «. — ► L 9 operator norm) is dominated by the geometric 

series which is convergent if -p- < 2 (— + l) S g , which is equivalent with q > q m . Therefore, 

ff= ft*,, 00^£ g PO, a>g m . (5.20) 

This finishes the proof of Theorem 13. Ill 
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6 Estimates for the half-wave operator 

Let (M, g) be a compact Riemann manifold of dimension n. Let P = \J — A + 1, where A is the Laplace 
operator. The principal symbol p(x, £) = g lJ (x)£i£j of P generates the Hamiltonian flow $ t : T*M — > T*M; 
this flow leaves invariant the cosphere bundle 

S*M = {(x,d)eT*M :p(x,0 = l}. 

The geodesies of unit speed on M are the curves t i— > 7r$t((x, £)), (a;, £) G S*M. Let 7r be the canonical 
projection T*M — > M. We say that the time t is non-conjugate if the bicharacteristics which start at the 
moment t = at any point x G M do not form caustics in time t, so that 7r$t : S*M — > M is of maximal 
rank: 

rankd(7r$ t | s , M )(£)=n-l. (6.1) 

Here £ is a point in the fiber S*M of the cosphere bundle at the point x. 

Assume that at t = T the map 7r$ t : S*M — > M is no longer of maximal rank at the point (x, £), where 

rankd(W t | s , M ) (£) <?i-l, t = T. (6.2) 

The integral kernel if* of the half-wave operator e ltp can be represented as a finite sum of oscillatory integrals 
of the form 

K t [x,y) = [ e« x *-« t <y'®at(v,t)dt, 

JR™ 

where at(y,£) is a classical symbol of order 0. (See |Sog93| , Section 4.) This representation is valid for 
(x, £, y) supported in a small open conic neighborhood of M x 1" x M and t in a small open neighborhood 
of T. We apply our results on LP — ► L 9 estimates (Theorem 13.5(1 to the half- wave operator e ltp with the 
integral kernel K t (x, y). 

Theorem 6.1 If for t < T the geodesic flow <I> t forms only caustics of the type A m i + i with m! < m, then 
for < t < T and for 1 < p < p < oo suc/i £/ia£ (p, q)^ ^ £ m i/ie LP — > L 9 estimates are caustics-insensitive. 
Precisely, 

e itPp-nS p -6 t . L p _^ L?) e 2tmj ( 6 .3) 

e ttp p -ns q -s p . L p ^ L?) ^ g jt e «g m . (6.4) 

For (p, qy G £ m , the estimates depend on the order of the caustic, which is given by n — | — m ^ 2 : 

elt p p - n i f -(« P +V(i/2+«)-(«,-y . ^ L 9 , (p, g )t G € m , q < p', (6.5) 
e itP p -nS g -(S p +S g )(l/2+K)-(S p -S 9 ) . L p ^ L q^ ( p;? )t £ Cmj g > (6. 6 ) 

T/ie regions 2t m , 2$ m , and £ m m (1/p, 1/q) -plane are defined in Definition \3.4\ (see also Figure^). 

We can use these results to investigate precisely the blow-up of the solution just before the formation 
of the caustics. At non-conjugate times t, the estimates on the half-wave operator e itp are given by the 
estimates (I6.5|) . (|6.6ll with k = 0. As t approaches the moment T when the geodesic flow starts forming 
caustics, these estimates blow up (and the estimates with nonzero n are to be used). As was shown in 
|Mag01| , if T > is such that t is non-conjugate for t € (T - e,T), for some e > 0, then the L q ' — ► L 9 
estimates on the half-wave operator e ltp may blow up as t — > T at most as 

\\e UP P- {n+1)/2 \\ Lq ^ Lq <C qM {T)\T-t\- K ^-^, 2<g<oo, (6.7) 

where K = 4. This is an a priori value; -ftT could be shown to be smaller when the geodesic flow forms some 
particular caustics. 
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Theorem 6.2 Let T > and suppose there exists e > such that t is non- conjugate for T — e < t < T . 
Assume that for T < t < T + e the geodesic flow <£"t forms only simple caustics of index at most m (e.g. 
caustics of the type A m+ i). Let 2 < q < oo, 1/q + 1/q' — 1. We have for T — e/2 < t < T: 

\W tP P' (n+1) HL^M)^(M)<C q .M{T)\T^t\- 5 ^ 2< g <oo, (6.8) 

We^P-^^-^h^M^LHM) < C q ,M(T), 2<q<^ (6.9) 

where S a = h — - and k = i t; . 

y 2 q 2 m+2 

Proof. We reduce the number of oscillatory variables in the representation of K f to 2, which is possible in 
an open neighborhood of simple caustics, and use the polar coordinates (A, a) £ R+ x § in the #-space. We 
exploit the fact that |<^ Q | > const |T — t\ if t is non-conjugate for T — e < t < T. (This bound is easy 
for stable caustics. For the generic situation, see |Mag 01 , Lemma 2.4.) We also use the Littlewood-Paley 
decomposition for K l (to interpolate L 1 — ► L°° estimates on K\ with L 2 — > L 2 estimates). The rest of the 
theorem is the same as the proof of the statement ()4.7Jl of Proposition l4.il Again, the optimal estimate with 
the factor | det<^ Q | -1 / 2 for the oscillatory integral is readily available since a is one-dimensional. □ 

The interpolation of the L p — > L q estimates which remain valid at the caustics (Theorem 16.1(1 and the 
asymptotics which describe the blow-up of the usual L p — > LP estimates (Theorem 16 .2(1 gives the complete 
description of the behavior of the blow-up of L p — > L q estimates just before the geodesic flow forms caustics. 

A Consistency of the definition of T> 

In this section we prove Lemma 12.71 



Lemma A.l (Lemma 12. 7|) Let A be a smooth closed conic Lagrangian submanifold of T*(X)\0. Let 
<p(x,8) € C°°{X x M. N ) be a smooth non- degenerate phase function which parametrizes A: 

A = {(x,d x cj)(x,e)) : d e (x,e)=0}. 

Let a = {cti] , 1 < i < N — 1, be local coordinates on the unit sphere S^ -1 . We use (A, a) £ M+ x S N_1 as 
local coordinates in R N . Then T> = dety (A -1 0^. Q J A ) ; 1 < i, j ' < N — 1, is a smooth function on A defined 
up to a nonzero factor: 

DeC°°(A)/^(A). 

We split the proof into two parts: In the first part, we will show that if we use the maximal number 
of oscillatory variables, then T> is defined up to a nonzero factor. In the second part, we show that T> is 
multiplied by a nonzero factor if we reduce the number of oscillatory variables. 

(z) Let us check that, up to a factor, T> does not depend on the chosen parametrization of A if we use 
the maximal number N = n of oscillatory variables. A can be parametrized (locally) by 9 with 9 £ t w , 
N = n. Assume there are two different phase functions 4>(x,9) and tp{x,d), 9 € 1^, £ R , and that 
both 9 and $ can be used as local coordinates on A. According to e.g. |Dui96| . there is a function g(x,9), 
homogeneous of degree 1 in 9, such that </>(x, 9) — ip(x,g(x, 9)). We rewrite </> and tp as 

4> = cj)(x,\,a), ip = ip(x,T,/3), (A.l) 

where A = r = and a, (3 are local coordinates on S N_1 . Then there is a smooth function (3(x, a) and 
a smooth function c(x, a) ^ such that 

<j>(x, A, a) = ip(x, c(x, a)A, (3{x, a)). (A. 2) 

To simplify the notations, we will assume the summation with respect to the repeating indices and will not 
write the subscripts of a, j3, ... at all, assuming that e.g. det A'^ a , stands for det^ d ai d aj A and A" aa i da da' 
stands for V, - d a ,d a .Adcti da-j. 
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We differentiate relation (|A.2|I twice with respect to a: 



where jf (a;, a) — d/3(x, a)/da (that is, jf ■ (x, a) — d/3j(x, a)/dai). Taking into account that ip' x \ = (ip is 
homogeneous of degree 1 in A), while ip x , ip'p, and ip'£p(x, A, /?) = A - ip'p(x, A, /3) vanish identically on A, we 
deduce that 

det<j)a a '(x,\,a) = (det jf (x, a)) 2 det tppp,(x, c(x, a)X, f3(x, a)), 

where det Jg(x, a) = det9 Q /3(a;, a) = det^ d ai (3j(x, a) ^ 0. 

(zi) Let us check that T> as an element of C°° (A) /C^ (A) is not affected by the reduction of oscil- 
latory variables. We consider the phase function <j>(x,\,a). Assume that the coordinates a split into 
a = (p, a) so that <$>" aa , is non-degenerate. Then there exists a smooth function S(x, p) such that the 
condition 4>' a (x, A, p, a) = is equivalent with a = E(x, p). The phase function 



ip(x, A, p) = <t>(x, A, p, S{x, p)) 



(A.3) 



parametrizes the same canonical relation as 4> does. We are going to prove that det $" aa i and det ip pp , differ 
by a nonzero factor (namely, det §" aa i ) • 

In what follows, we drop off the dependence on x and A. Differentiating (|A.3(I with respect to p, we get 



V>» = </>>, S{p)) + ^(p, E{p))r p {p), (A.4) 

where J a p {p) = dS(p)/dp. 

KM = € P '(P> S (P)) + C(A ApWM + <#„(p, m)Jp(p) (A.5) 

+KAp, s{p))j;{p)j;;{ P ) + <^(p, £(p))j; p ,(p). 

The last term in the right-hand side of (|A.5(1 vanishes identically on the canonical relation (where <p' a = 
((/) p , ftcr) = 0). Using the identity 

= d p (^(p, 2{p))) = <t>%(p, S{p)) + C,< (p, S{p))j;\p), 

we can expr ess J p (p) = — 0" , (p, E{p))<p aa (p, S{p)), where <T<T (p, cr) denotes the matrix inverse to <$>" aa , (p, a). 
We rewrite HA.5[) as 

- z(p)) - €M Ap)W°' (p, z{p))CAp, z(p)). 

To compute the determinant of l|A.6j) . we use the identity 



(A.6) 



A B 
C D 



det {A - BD- l C) det D = det 
where A and D are square matrices and det D ^ 0, which follows from the matrix identity 



(A.7) 



A-BD- l C 







' I 


-BZ)- 1 " 




' A 


B ' 




I " 





D 







7 




C 


D 




—D~ X C I 



Identity (|A.7|) allows to write the determinant of (|A.6(1 in the form of the desired relation: 

det J," det Aoo> = det < Q , , where det <^ CT , ^ 0. 



(A. 



This finishes the proof of Lemma \l. 71 
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B f) 1 — > L°° estimates on (A, cr)-pieces 

In this section we prove the following Lemma, needed for sharp f) 1 — > L q and i p — > BMO estimates. 

Lemma B.l Let ^ G I^iX, Y, C) fee associated to a canonical relation such that C -> X is a submersion, 
and let C /iaue only caustics of the type A rn+ i with m = 1 or 2. Then, for any atom (Zq supported in the 
cube Q wi/i side r, we /iai>e 

||A-"-( n+JV )/ 2 ffA,aaQ|U- W < const^/'-minlAr^Ar)- 1 ). (B.l) 

For simplicity, we consider $ £ I^{X, Y, C) with // = — (TV + n) /2. This implies that a(x, 8, y) £ S d with 
d = — (n + N)/2 — (N — n)/2 = —N. Let oq be an atom supported in the cube Q with side r (following 
SSS91], we may assume that r < 1). We want to show that for any x, 

\3\,*a Q (x)\ < const cr 1 / m min(Ar,(Ar)- 1 ). (B.2) 

We will decompose and bound the pieces 5a,o- following the discussion on pages 238-241 in SSS91 . For 
a particular A, we introduce unit vectors A , with 1 < v < N(\~ 1 / 2 ) m X~ ~ , equidistributed on the unit 
sphere in the 0-space M. N , so that \9% — 9\ | > const A~2 for v ^ v' . We introduce a corresponding partition 
of unity, 

1= E WW), 

where the functions %1>X are homogeneous of degree and supported in the spherical angles f2 A with the span 
~ A~ 1/2 , centered at 9K: 



Vx(0)^O only if 



1*1 



< const A 2 . 



We assume that \d^(9)\ < const Att |0|-M . 

We denote the integral kernels of 3\, a by K\ :Cr (x, 9, y) and K\ t<r (x, 9, y). We introduce 3a a by 



$\, (T u{x) = J K v X(J (x,9 1 y)u(y)d9dy, 

where K^x, 6, y) = ^ (9)K Xt<7 (x, 9, y). 

From now on, we assume that x £ X is fixed. We need to introduce the "exceptional set" associated to 
x. According to |CC03j . the assumption that C — > X is a submersion allows one to choose the phase <ft in 
the form 

<Kx,6,V) = (G(x,9,y")-y',6), 
where y = (y', y") £ R N x R" _JV = Y are certain local coordinates. For given A and v, we define TV A _i C Y 

by 

= iv ■■ \( G ( X > d x>v") - y'> < A ^ W*, d hv") - v'\ < A ~n, (b.3) 

with [R, v x A -i I < const A -1 • A ~ . We set xtz v _ x ( x ) to be the characteristic function of 1Z V X A _ t . 
Let aQ be an atom supported in the cube Q with side r: 

|Q| = r", ||a Q || L »<r- n , ||a Q || £1 < 1, [a Q = 0. 



Q 



We consider 



3\,aa Q (x) = E/ Xk» x _ 1 {y)KZ a {x,9,y)a Q (y)d9dy 

+ J2[( 1 -^:, x -M^ K l^ d 'y) a Q^ d6d y- ( R4 ) 
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We need to know the absolute value of this expression. 

(i) The absolute value of the first term in the right-hand side is bounded by 

/ Xv.^ x _ 1 (y)\K^ CT (x,e,y)a Q (y)\ d6dy 

< C\- N J2 I XTz :iX _M^( 9 W\/X)0(V(x'',e,yya)\aQ(y)\dedy. (B.5) 

o In 1|B.5() . we have already applied the bound CX~ N on the symbol a{x,8,y) E S~ N at \8\ ~ A. 
o Summation in v converges since ^„ ip\(0) = 1. 

o If dgT) ^ (m = 1), then the integration in 8 contributes const aX N , where a appears due to the 
support properties of (iiV/u). 

More generally, assume that at a point p S C there is a simple caustic of the type A m+ i with m = 1 or 2. 
Then there is a vector field V = ajd 6j , V G C°°(T(T(C))), such that V m V\ p ^ 0. We define 9 = f e R N , 
so that the region of integration in 8 is bounded uniformly in A. Note that d9 = X N d<d. We can choose 
the coordinates so that dg l N T> ^ 0, in an open neighborhood of p. The expression is homogeneous of 

degree zero in A, so that \dQ N V\ > const > uniformly in A, a. Therefore, 



/ (3(V/a)de N < const cr™ 



The integration in ©i, . . . , &n-i converges since the support of IB. 5(1 in O = j is bounded (uniformly in A, 
a). We conclude that the integration in 8 contributes const \ N <j™. 

o Finally, due to the bound HaQllioo^y) < |Q| _1 together with the support properties of aQ and 
Xw _! the integral in y contributes the factor min(l, (Ar) _1 ). 

Taking the product of all of the above factors, we obtain const <r™ min(l, (Ar) _1 ). 
(ii) For the absolute value of the second term in the right-hand side of l|B.4(l we have: 
o In each v-teim, we can integrate by parts as in |SSS91j (we need the assumption a > A~2 to obtain 
an analogue of the inequalities (3.19) in SSSiH]; the argument is the same as theirs), getting the factor 

(1 + A 2 | (G(x, 6 v x ,y") - y', 9» x ) | 2 + X\G(x, ^, y") - y'\ 2 y M , for any M e N. 



The integral of the product of this expression with csq (y) with respect to y contributes the same factor 
min(l, (Ar) _1 ) as above. The rest of the analysis is the same as for the first term in the right-hand side of 
(lB~4ll . 

We conclude that |3A,<7<iq(2;)| < const am min(l, (Ar) _1 ). 

The bound in the case Xr < 1 follows from |CC03j . Let us recall the argument. We fix some point y S Q. 
Since J aQ(y)dy = 0, we can write 

Sxo-aQ^) = f [K\ <a (x,6,y) - Kx,a{x,9,y)]aQ(y)d6dy 
i 



dtdt^J K x ,a(x,8,y+ {y - y)t)a Q {y)dd dy 
= Xr J (J dt^^X~ 1 d y K x . a (x,8,y+(y-y)t) >l ja Q (y)d8dy. (B.6) 

The expression in the curly brackets can be treated as an integral kernel of another Fourier integral operator 
of the same order \i associated to C, and therefore 



l/m 



\-&\,aO<Q(%)\ < Xr const a 

Let us mention that in i|B.6() | | < const and that the increase in the order of the symbol due to the 
derivative d y is compensated by A -1 . When the derivative d y acts on f3(D(x, 8, y)/u) (which is hidden inside 
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K\ t(T ), the contribution is bounded by const c -1 and is also compensated by A -1 . The integration in t is 
irrelevant. 

This completes the proof of Lemma IB. II 
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